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Abstract 

The recently observed color suppressed B° — > D ir°, D rj^'\ DfK~ and D°K° decay modes 
all have rates larger than expected, hinting at the presence of final state interactions. We study 
rescattering effects in B — ► DP, DK and DK modes in the quasi-elastic approach, which is 
extended to accommodate D°n' without using U(3) symmetry. The D°K modes are of interest 
in the determination of the unitarity angle ^3/7. The updated DP data are used to extract the 
effective Wilson coefficients af ~ 0.92, a| — 0.22, three strong phases 5 ~ 62°, 9 ~ 24°, a ~ 127°, 
and the mixing angle r ~ 2°. The values of 5 and 9 are close to our previous results. The smallness 
of r implies small mixing of D°rji with other modes. Predictions for D°K~, D + K~ and D°K° 
agree with data. The framework applies to B — > DK, and rates for D°K~, D~K°, D~7r°, D~n 
and Djrf modes are predicted. From B~ — ► D°K~ and D°K~ rates, we find rs = 0.09 ± 0.02. 

PACS numbers: 11.30.Hv, 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



The color-suppressed decays B — > D^*> n 

HQ 

and D°r], D°uj |l[ were observed for 
the first time in 2001. Recently, improved measurements of 5° — > D^°(7r°,rj,uj) 0, |^ 
and the first observation of the D°r]' mode have been reported by BaBar [sj] and confirmed 
by Belle 5 . Other color suppressed modes, such as D S K~ and D°K°, have also been 
observed 6, Q|. All these modes have branching ratios that are significantly larger than 
earlier theoretical expectations based on naive factorization, indicatin g the p resence of non- 
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vanishing strong phases, which has attracted much attention 

\v\ . . Shortly after the first observation of the color suppressed modes became known, we 
proposed ^| a quasi-elastic final state rescattering (FSI) picture, where the enhancement 
of color suppressed D°h° modes can be understood as rescattering from the color allowed 
D + 7i~ final state. The framework is applicable to B — > DK, D K decays. 

The color-allowed B~ — > D°K~ and color-suppressed D°K~ decays are of interest for 
the determination of the unitary phase angle 03 (7) = argV^, where V is the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix. To be more specific, the amplitude ratio 
r B and the strong phase difference 5 B for D°K~ and D°K~ decay modes, which are governed 
by different CKM matrices as depicted in Fig. Q are defined as 
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The weak phase 03 is removed from A(B — > D°K ) in defining 5b- The r B and 5 B param- 
eters are common to the 03 determination methods of Gronau-London-Wyler (GLW) [19], 
Atwood-Dunietz-Soni (ADS) Q and U DK Dalitz plot" 0, Q> where one exploits the 
interference effects of B~ — > D°K~ — > fcpK~ and B~ — > D°K~ — > fcpK~ amplitudes. 
Note that the r# parameter, which governs the strength of interference, is both color and 
CKM suppressed, hence hard to measure directly. 

Through the DK Dalitz plot method, the BaBar and Belle experiments a 
7 = 70° ± 44° ± 10° ± 10° and 3 = 64° ± 19° ± 13° ± 11°, respectively 0, 
the last error comes from modelling of D decay resonances across the Dalitz plot for, e.g. 
D° — > Ks^ + ^ ■ Although similar results on 3 are obtained, the corresponding values 
are quite different for BaBar and Belle. Belle reports = 0.21 ± 0.08 ± 0.03 ± 0.04, while 
BaBar gives r B < 0.19 at 90% confidence level. Note that an average r B = 0.10 ± 0.04 
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FIG. 1: Color-allowed and color-suppressed B 
B~ -> D°K- decay. 



is found by the UTj it group, by combining analyses using all three methods |25j]. As the 
strength of interference is governed by the size of rs, the larger error in the 7 value of BaBar 
reflects the smallness of their r«. Given the present experimental situation that Belle and 
BaBar have quite different tb values and the critical role it takes in 0^3/7 extraction, it is 
important to give a theoretical or phenomenological prediction of rg and 5b- Not much 
work has so far been done. 1 

The enhancement in rates of the color-suppressed DP modes could imply [15] a larger r r . 



It is thus of interest to study DP and D K modes together. In fact, it was noted in Ref. 
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1 In the preparation of this paper, we note that a calculation in the pQCD approach has been reported [2^ . 
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that the quasi-elastic approach used in the DP system can be applied to the D K system. 
Rescattering parameters are basically non-perturbative and can only be fitted from data. 
We shall use the B — > Dn, Dr)('\ D S K decay rates to extract the rescattering parameters, 
which are then used to predict DK and D K rates. In this way, we are able to deduce what 
value of tb is preferred in the quasi-elastic rescattering scenario. 

Our previous analysis on FSI effects in D°h° modes was based on early data. Recent 
experimental updates show some variations. For example, the D°ir rate has dropped while 
the .Derate is larger. Furthermore, the D°7]' mode has finally been measured. In our earlier 
study (ill, because of the absence of the D°r]' mode, we ignored it and approximated D°r] 
by D°i]s, by argument of the Ua(1) anomaly and small singlet-octet (or 77-7/) mixing. The 
same approach was applied to the study of the charmless case j^. Given the long standing 
problem of the B — > rjK rate, it is of interest to clarify the r\\ issue in B decays |2£j. With 
the emerging new data, it is time to update the quasi-elastic rescattering approach, and 
verify the approximations made. 

In Sec. II we extend the quasi-elastic rescattering framework to include D°r]'. Because of 
the Ua(1) anomaly, we use SU(3) rather than U(3) symmetry. As a consequence, we need 
three phases and one mixing angle as rescattering parameters. These same parameters also 
enter the rescattering in the DK and D K systems. In Sec. Ill we carry out a numerical 
study. The effective Wilson coefficients and rescattering parameters are obtained by using 
current B — > DP data. The B — > DK rates are then predicted and compared with current 
data. We then proceed to study the B — > D K system and make predictions for r B and 8b- 
The conclusion is then offered in Sec. IV. 

II. FINAL STATE RESCATTERING FRAMEWORK 
A. Quasi-elastic Rescattering 

Let if\y be the weak decay Hamiltonian. In the absence of weak phase (or if they are 
factored out), if\y is time- reversal invariant. By using time reversal invariance of if\y and 
the optical theorem, we have (see, for example, [ill I29I]) 

2Im(i;out|# w |£) = ^ 7}* (j; out|F w |£>, (2) 
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where T is the T-matrix of strong scattering, and the phase convention T |in) = |out) under 
time-reversal operation is used. This is the master formula of FSI in B decay. In particular, 
for B decay to two body final state with momentum (pi,P2), we have, 

" 2i1mA(p B - PlP2 ) = E (fl/ jS^) (27F)V ( Pl+P2 " 2> 

j ' \h — 1 / \ k — 1 

xM*(p lP2 - {q k })A(p B - {g fc }), (3) 

where the optical theorem is used to all orders of the strong interaction, but only to first 
order of the weak interaction. Eq. (j3J) relates the imaginary part of the two body decay 
amplitude to the sum over all possible B decay final states {qk}, followed by {q k } — > pip 2 
rescattering. The solution to the above equation is 

A = S l ' 2 A\ (4) 

where A is real, and S = 1 + %T . The weak decay picks up strong scattering phases (^(J. 

It has been pointed out that j3l| elastic rescattering effects may not be greatly suppressed 
at mg scale, while inelastic rescattering contributions may be important. But we would 
clearly lose control if the full structure shown in Eq. (J5J) is employed. Even if all possible B 
decay rates can be measured, it would be impossible to know the phases of each amplitude. 
Furthermore, we know very little about the strong rescattering amplitudes. However, in 
Eq. Q the subset of two body final states that may be reached via elastic rescatterings 
stand out compared to inelastic channels. It has been shown from duality arguments [32], 



as well as a statistical approach 29|, that inelastic FSI amplitudes tend to cancel each other 
and lead to small FSI phases. We shall therefore separate {qk} into two body elastic channels 
plus the rest, and concentrate on the contribution of the former. 

We will consider DP final states, where P stands for a pseudoscalar meson. We consider 
D and P within SU(3) multiplets, for example the D anti-triplet of D°, D + and Df, the n 
octet that contains pions, kaons and the r\ (rjg component) meson, as well as the rji (mixing 
of physical rj and rj'). Thus, we call this quasi-elastic rescattering. 

B. Rescattering Formalism including D°rji 



In Ref. we treated -DII — > DU rescattering and considered only the rj^, since at that 
time B° -> D° rj' was not yet reported. Here we wish to extend the formalism to include 



D°r]i, so we can treat the physical D°r) and D°r]' final states. The quasi-elastic strong 
rescattering for C — +1, S — final states can be written as 

/ A D+W - \ 



where S 1 ' 2 = (I + iT) 1 ' 2 
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and T has the same structure as T by SU(3) symmetry, but with replaced by r[. In 
addition, we have A D o n - = (1 + ir' Q + ir' e )A° D07r _. The rf s are discussed below. 

Eq. (jUJ) is obtained by separating the DU — > DU scattering amplitude M into three 
independent components, M a e , and defining 
d 3 qi d 3 q 2 



(2vr) 5 (pi +p 2 -qi- q 2 ) Mi(p x p 2 -> q x q 2 ). 



(7) 



(2tt) 3 2.Ei (2tt) 3 2£ 2 

Note that the Mo, M a and M e amplitudes correspond respectively to the three independent 
SU(3) combinations Tr(D ill Z) out )Tr(ri ili ri out ), Tr(A n n in n outJ D out ) and ^(DiJI^IL^D^) of 
D in Il in — > _D out n out scattering. For DU <-> D°r]i and D°rji <-> D ?]! scattering, we denote the 
corresponding integrals as fj and fj, respectively. Had U(3) rather than SU(3) symmetry 
held in the light pseudoscalar sector, fj and fj would have been identified with rj. 

We give a pictorial representation of r e , r a , r in Fig. El which can be seen as corresponding 
to charge exchange, annihilation, and flavor singlet exchange, respectively. The coefficients r« 
in Eq. (JBJ) can be reproduced easily using this pictorial approach by matching the flavor wave 
function coefficients. For example, we have (r a — r e )/\^2 for D + n~ — > D n° rescattering. 
Exchange rescattering (r e ), the first diagram of Fig. projects out the dd component of 7r° 
on the right hand side. This give a —l/y/2 factor from the it wave function. Similarly, the 
second diagram of Fig. [21 projects out the uu component of 7r°, hence gives r a /y/2. These 






FIG. 2: Pictorial representation (from top to bottom) of charge exchange r e (r e ), annihilation r a 
{f a ) and singlet exchange r$ (fo) for DP (re)scattering. 

diagrams also provide further information. For example, it is easy to see from the second 
diagram that annihilation rescattering (r a ) is responsible for D + n~ — > D+K~ , since there 
is no s quark before rescattering. The MjS in Eq. (J7J) can be understood analogously. For 
example, the M e operator (Di n )i(n out ) l k (Ili n )j(D ont y corresponds to 2 the D(ccji)Il(qjq k ) — > 
D(cqi)~n.(q l qk) process with the exchange of ith and jth antiquark, hence it is the exchange 
rescattering operator. 

For Cabibbo suppressed B — > D + K~ and D°K° modes, we have 

A D+K -\ = fl + ir' ir' e W A° D+K 
A d oko ) \ ir' e l + ir' J \A° D0W 

2 Superscripts and subscripts are assigned according to the field convention with q l and ijj as quark and 
antiquark fields, respectively. 
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FIG. 3: Pictorial representation of charge exchange r e in D+K~ D°K° (re) scattering. 



which can be easily obtained by using the pictorial approach shown in Fig. El It is clear 
that "annihilation" rescattering is impossible in this case. For D°K~, we have A d ok- = 
(l + ir' + ir' e )A DOK ^. 

For the B~ — ► D K decays, we have 
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(10) 



r + 3 

as one can easily verify using the pictorial approach shown in Fig. |3J The zeros are a 
consequence of assuming isospin symmetry. It is important to note that, due to charge 
conjugation invariance and SU(3) symmetry of the strong interactions, these rf , and 
r^P coefficients are identical to those in Eqs. © and (JBJ). 
By solving S^S = 1, we obtain [h| 
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FIG. 4: Pictorial representation of charge exchange, annihilation and singlet exchange 
(re)scatterings involving the D°K~ state. 



(! + *■„) = \{l + e 2 % 

%r a = h3Un-2e 2iS -l), 

o 

/- 3 

i[r a + r e ) = ^Wi 2 , 
^(ro + !: ^) = Wjb-1, (11) 



where 

„ / cos r sin r \ / e 2ie \ / cos r — sin r , 
U = U T =\ | . (12) 

— sin r cos r / y e / y sin r cos r 
where we have set the overall phase factor (l + iro + ir e ) in 5 to unity. This phase convention 



is equivalent to choosing the A D o n - amplitude to be real. 

We stress that the above solution satisfies S^S = 1 in all three cases of DP (P now stands 
for strangeness pseudoscalar) , DK and DK. The r-, f[ and f[ in S 1 ^ 2 can be obtained 
by using the above formulas with phases (5, 8, a) reduced by half. We need three phases 
and one mixing angle to specify FSI effects since one does not have nonet symmetry in the 
light pseudoscalar sector. An extra phase as well as mixing angle arise from including rjx in 
our analysis. We will resort to data to see how far the DP system differs from the nonet 
symmetric case. At the same time, we will use r-, f- and to predict DK and D K rates, 
and compare with data whenever possible. 



C. SU(3) Decomposition 

It is instructive to see the phases and angle given in Eq. (|11|) in light of SU(3) decom- 
position. Let us consider the DP case first. D is an anti-triplet (D(3)), while P can be 
reduced to an octet [11(8)] and a singlet (r/i). The D(3) (g) 11(8) can be reduced into a 3, a 6 
and a 15 (see, for example, j^), while D(3)r]i is another anti-triplet. Denoting the latter 
as 3', it can mix with the 3 from DU via a 2 x 2 symmetric (from time reversal invariance) 
unitary matrix U, which appears already in Eqs. (jllj) and (fT2*|) . The invariance of the strong 
interaction under SU(3) transformation gives 

S= |T5)(15| +e M |6)(6| + (|3> \3')) -U ■ \ _ ). (13) 

V < 3 'i / 

It is now clear that, with the choice of vanishing phase in Sj^j^, 25 is the phase of and U 
is the mixing matrix in the anti-triplet sector. Note that in the master formula Eq. (j3J) one 
should use S 1 ^ 2 . This can be easily obtained by reducing all phases in the right-hand-side 
of the above equation by half. 

Three remarks are in order. It is important to emphasize that, by charge conjugation 
invariance of the strong interaction, the above ^-matrix can also be applied to the DP case 
with 15, 6 and 3^ replaced by 15, 6 and 3^, respectively. Second, the D°r] 8 and D ^ are 
not physical final states. The physical t], r\ mesons are defined through 
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(14) 



with the mixing ang letf = -15.4° Q In the previous analysis [llj , D rj was approximated 
as D°7]s, while D°r]i was ignored. This corresponds to vanishing mixing angle r = in U, 
and dropping the 3'. 

Before we end this section, let us specify A in the above formulas. We use factorization 
amplitudes for A , which are 0|, for each physical final state 3 



4o.- = V cb V: d (T f + C f ), A f D+n _ = V cb V* d (T f + E 
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where the super- and subscripts / indicate factorization amplitude, and 
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(16) 

F BD ( BF "> i s the B — > D(P) transition form factor and Fq^ dp is the vacuum to DP (time- 
like) form factor. Some SU(3) breaking effects are included in the decay constants and 
various form factors. In subsequent numerical study we take Ef = af = 0. In this case, it 
is seen that in the B DP decays, the af and af parameters accompany F BD (0) and 



3 To match the phase convention T |in) = |out) (after factoring out the CKM factor) we absorb the factor 
i of these factorization amplitudes into the B meson state |l8j. 
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TABLE I: Summary of experimental results for B° — > DP modes 



QBQBBBQ 



B (xlO" 4 ) 


CLEO 


BaBar 


Belle 


Average 


5° -> D°7] 


2.74lg;| ± 0.55 


2.9 ±0.2 ±0.3 
2.5 ±0.2 ±0.3 


2.31 ±0.12 ±0.23 
1.83 ±0.15 ±0.27 


2.53 ±0.20 
2.11 ±0.33 a 


B° -> D+K~ 
W -> D r/ 
5° -> D + K~ 




0.32 ±0.12 ±0.08 
1.7 ±0.4 ±0.2 


0.45^2 + 0.11 
1.14±0.20±g;^ 
2.04 ±0.50 ±0.27 


0.38 ±0.13 
1.26±0.23 6 
2.0 ±0.6 


W -> D°l< 






0.50±g;i| ± 0.06 


0.50 ±0.14 



"The error is scaled by a factor of 5 = 1.4. 
6 The error is scaled by a factor of S = 1.1. 



TABLE II: Form factors in covariant light-front models |4J|. For — > r/ form factors the mixing 
angle and Clebsch-Gordan coefficients are included. 



in 



D,D S 



7 Br, 




FP K {m 



2 ) 
D,D.) 



0.28 
0.15 
0.43 



rpBD 
^0 



0.67 
0.13 



Fq p {m 2 D ) 1 respectively. It should be noted that in the quark diagram approach [or the SU(3) 
approach] the full amplitudes (i.e. without the super- and subscripts /) are expressed as in 
Eq. (JT3J), and one treats T (tree), C (color -suppre ssed), E (exchange) and A (annihilation) 



as complex topological amplitudes |35j, |3£ , 



37, 



38, 



III. RESULTS AND DISCUSSION 



In our numerical study, masses and lifetimes are taken from the Particle Data Group 
(PDG) 40]. We use the color suppressed branching ratios as stated in Table UJ For other 
modes, such as B° — > D + n~ and B~ — > D°ti~ , D°K~ decays, we use PDG values J^. We 
fix V ud = 0.9738, V us = 0.2200, V cb = 0.0413, V cs = 0.996, \V ub \ = 3.67 x 10" 3 , and use 



the decay constants f n = 131 MeV, f K = 156 MeV Q and f D = 200 (230) MeV. Form 



factors are taken from the covariant light-front quark model calculation [4l| . where we list 
the relevant values in Table HU 
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To describe the processes with rescattering from factorization amplitudes, we have six 
parameters: the two effective Wilson coefficients af and af, and the three rescattering 
phases 5, 9 and a, and one mixing angle r in S 1 ^ 2 . In subsection A, these parameters are 
fitted with C = 1, S = DP data, i.e. the B -> £>+tt-, D°tt-, D°7r°, D\ D°rf and D+K" 
rates. We then use the extracted parameters to predict (C — 1, S — — 1) B — > D°K~, 
D + K~ and D°K° rates and compare with measurement. Predictions for the (C = — 1, 
S — — 1) B ^ D K modes, and the value for r#, are given in subsection B. 

A. FSI Effects onB-t DP and L>K Rates 

Taking the B -> D+tt", D ^, D°7r°, Z^r/, D°r/' and rates as input, we fit for af 2 

and the FSI phases and mixing. The fitted parameters are given in Table IHIl where errors 
are propagated from the experimental errors. The factorization rates and the predicted 
DK rates are compared to experimental results in Table IIVI The factorization rates are 
determined by using the af 2 from the fit, but setting all FSI phases and mixing angle to 
zero. Unitarity is then implied automatically, i.e. sum of rates within coupled modes are 
unchanged by FSI. 

Table ITVl illustrates the effect of FSI. B° — > D°h° (where h° = ir°, 77, 77') rates are fed 

TABLE III: Fitted parameters of the SU(3) FSI picture, using B D + n~, D°ir-, D ir°, D°rj, 
D°7]' and D^K~ decay rates (see Table ITV)l as input. There is a two fold ambiguity (the overall 
sign of the phases) in the solutions. The SU(3) phases and mixing are reexpressed in terms of the 
rescattering parameters r' t , f[, f\. 



parameter 


solution 


parameter 


solution 


af 


n Q2+ - 04 


af 


22+ 012 


5 


±(62.411;*)° 


e 


±(23.8^;°,)° 


a 


±(127j^)° 


T 


(1.7tS 8 )° 


l + ir' 


(0.73 ±0.04) ± (0.44±0.02)i 


K 


(0.27 ±0.04) =F (0.44±0.02)i 


ir 'a 


(0.10 ±0.02) T (0.07±0.02)i 


i(r' a + K) 


(-0.051°;^) ± (0.01±g:i;)i 


1 + ir' Q + i T -^^ L 


(-0.61^) ± (0.79±?^)i 


iK + <) 


(0.37 ±0.06) =F (0.51 ±0.03)? 


1 + ir' + i T -^j^ 


(0.86 ±0.02) ± (0.27ig;^)i 
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TABLE IV: The branching ratios of various B — ► DP and DK modes in 10~ 4 units. The second 
and third columns compare experiment with factorization model, and the last column gives the 
FSI results. The factorization results are obtained by using the same set of parameters but with 
FSI phases set to zero. 



Mode 


K exp (io~ 4 l 

a {±v ) 


K fac (4 0~ 4n I 

ID \1U j 




D°7T- 


49.8 ± 2.9 


input 


input 


D+TT- 


27.6 ±2.5 


33.0+1^ 


input 


D°ir° 


2.53 ±0.20 


U - O1 -0.34 


input 


D+K- 


0.38 ±0.13 





input 


D°7] 


2.11 ±0.33 


0.29^ 


input 


D°rj' 


1.26 ±0.26 


o.i8i8:S 


input 


D°K- 


3.7 ±0.6 


o ni+0.37 


o ni+0.37 
,3 - yl -0.32 


D+K- 


2.0 ±0.6 


9 OO+0.21 

^• oo -0.31 


1 78+ - 20 

L - '°-0.17 


D°K° 


0.50 ±0.14 


12+ 017 


n 7 o+0.08 
u - '°-0.10 



mostly by B° — > D + 7r _ . Since these rescattering parameters are extracted from CP even 
measurements, the overall sign of phases is undetermined. 

The FSI contributions for the D^K~ , D°tt° and -D°?78 rates from D + tt~ rescattering are 
governed by r' a , r' a — r' e and r' a +r' e , respectively. The strength of B(B° — > DfK~) ~ 4 x 10~ 5 
implies that r' a cannot be too small, i.e. \r' a | ~ y/B(D+K-)/B(D+Tr-) ~ 0.12. On the other 
hand the FSI enhances D°ir° rate from 0.5 x 10~ 4 to 2.5 x 10~ 4 . Comparing these two 
modes, we have \r' e — r' a \ > \r' a \. Consequently, through the analysis of the two above modes, 
the size of FSI contribution to -D°?78 is roughly determined. In fact, the FSI contribution 
alone already gives B(D°r]g) ~ 2 x 10~ 4 , and after interference with the short distance 
contribution, one gets B(D r] 8 ) ~ 3 x 10~ 4 . To fit the D°r] data, the D°r] 8 amplitude 
interferes destructively, through a > 90°, with the D°r]i amplitude, while at the same time 
the D°r]' amplitude gets enhanced through constructive interference. Although the D°i]i 
amplitude is small and does not be enhanced in FSI (due to the smallness of f' ea ), it still 
affects D°T] and D°r]' rates through interferences. 

Utilizing the SU(3) framework, we can predict the results for C = 1, S = — 1 DK 
modes. The prediction for B — > D°K~, D + K~ and D°K° rates are given in Table ITV| 
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where the experimental data are also listed. These rates were not used in the fit for at and 
FSI parameters. We see that the D°K~, D + K~ rates are in good agreement with data. 
The color-suppressed D°K° rate is a bit larger than data, but still in reasonable agreement. 

From Table ITTT1 we observe that 1 + ir' ~ 0.85 e ±l31 ° is almost perpendicular to ir' e ~ 
0.52 e^* 58 ". We do not know the reason for this orthogonality, but this implies that the 
FSI amplitude from D + K~ — > D°K° rescattering is almost perpendicular to the one from 
D°K° — ► D°K° rescattering [c.f. Eq. (jHJ)]. Consequently, we have 

B(D°K°) ~ |1 + ir' \ 2 B iac (D°K°) + \ir' e \ 2 B^ c (D + K~), (17) 

which gives a very good approximation of the result shown in Table IIV1 As we shall see, a 
relation similar to Eq. (J 17)) holds for the B~ — > D°K~ case. From Eq. ()17|) we see that a 
~ 15% reduction of a\ s r' e from its central value can reproduce the current B(D°K°) central 
value. In fact, a smaller B(D°K°) ~ 0.5 x 10~ 4 was predicted [11] with a smaller 5 extracted 4 
from earlier D°h° data, so the measurements probably have yet to settle. Note that the rate 
of the color- allowed D°K~ mode is not affected by the quasi-elastic rescattering, just like 
D°n-. 



B. B — > DK Rates and Prediction of rs 



In this subsection, the B — > D K rates and Tb, &b are predicted and compared with data. 



Table 
method 



V] gives the current experimental results on re, 5b and ^3/7 from 5 the DK Dalitz 
. Our predictions for B — > D°K~, D~K°, D~tt°, D~rj and Djrj' decay 



TABLE V: Summary of experimental results for rs, $b and 4>3 in the DK Dalitz method 231. \24 ] 
For the phase convention adopted see footnote 03 



Belle 



BaBar 



tb 



^3(7) 



0.21 ±0.08 ±0.03 ±0.04 
-23° ± 19° ± 11° ± 21° 
64° ± 19° ± 13° ± 11° 



< 0.19 (90% CL) 
-66° ± 41° ± 8° ± 10° 
70° ± 44° ± 10° ± 10° 



4 In place of 8, a different notation 8' was used in 

5 Note that we use the phase convention CP \D°) = — \D°) and, consequently, our 8b is related to those in 
0,0 



by a 8b — tt transformation. 
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TABLE VI: Predictions for B~ —> D K rates. Experimental limits 40] are shown in the second 

column. The third and fourth columns are factorization and FSI results, respectively, using the 
same parameters as Table ITTT1 

Mode £ cxp (1(T 5 ) £ fac (1(T 5 ) B FS1 (10" 5 ) 

D°K- 0-17lo:ii °- 28 iai5 

D-K° 0.05±g;g| 

D 7 n° < 20 0.77^;?o 7 0.59^ 

D 7v < 50 0.46±g:8S 0.17±85§ 

Djr]' 0.30 ± 0.03 0.58±g;^| 



rates are given in Table IVII The sum of rates is unchanged in the presence of FSI. Note 
that the first two modes have I = 0, 1 components, the third mode is purely 1 = 1 while 
the last two modes are purely 1 = 0. Rescattering between D~n° and D~ry^ is forbidden 
by isospin. Consequently, the D~i] and D~rj' rates do not receive any contribution form 
the D~7T° mode. These two modes also do not rescatter much from D~K° and D°K~, 
as the rescattering parameters are either suppressed by Clebsch-Gordan coefficients, or by 
the smallness of f' e + f' a (c.f. Table HLT|) . Thus, these two modes basically rescatter among 
themselves (as one may check that the sum of their rates are roughly conserved under FSI). 

The D~rj and Djrf rates are reduced and enhanced, respectively, through FSI between 
themselves. This is due to the destructive and constructive interference effects of D~r] 8 and 
D~r]\ in the a > 90° required from the D°r}^ data. Note that within error we can 

also have a < 90°, so B(Djrj) > B(Djrf) is not ruled out. 

For the first three modes, the dominant source of rescattering is D~tt°. It feeds D°K~ 
and D~K° through r' e /\/2 and —r' e /^/2, respectively. Note that these S 1 ^ 2 matrix elements 
are similar to the D + K~ — > D°K° rescattering matrix element except for the l/y/2 factor. 
Furthermore, we see from Table IIHI that 1 + ir' Q + ir' a ~ 0.97e ±?32 ° is perpendicular to 
ir' e ~ 0.52 e Tl58 °. Consequently, the FSI amplitude from D~n — > D°K~ rescattering is 
orthogonal to the one from D°K~ — > D°K~ rescattering [c.f. Eq. (J5J)], and the rate of 
D°K~ is roughly given by 

2 



B(WK~) ~ B iac (D°K-) + 



K 
V2 



B iac (Dj7T°) ~ B tac (D°K~) + 0.13 B^ c (D;n°), (18) 
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which is analogous to Eq. (fTTjl for the D°K° case. Since the D°K° and D°K~~ modes have 
similar relations and the FSI contributions are governed by the same parameter ir' e , the 
D°K° mode may provide some estimation of the FSI effect in the D°K~ case. For example, 
if the current central value of 5 x 10~ 5 for D°K° holds, we would need roughly a 15% 
reduction in af^r' e as discussed in the previous section. This corresponds to a 10% reduction 
in the predicted D°K~ rate or, equivalently, a 5% reduction in |Aqo#-|. Such a variation is 
within the error shown in Table IVI1 

The B~ — > D~K° decay is a pure annihilation decay mode. Its decay rate can give us 
an idea of the size of the annihilation amplitude. We set the short distance annihilation 
contribution to zero. Its rate then comes mainly from D~7r° rescattering, which interferes 
destructively with the rescattering contribution from D°K~. In terms of topological am- 
plitudes we have (a/^^ = 0.21 ± 0.08 and arg(a/t)-^^ = ±(104.87Jia2)°> where a and t 
are the full annihilation and tree amplitudes, respectively [see Eq. ffTKjl and subsequent dis- 
cussion]. The (a/t)j)K ratio is numerically close to a| ff /a^ ff . The long distance annihilation 
amplitude cannot be neglected. 

The A(B~ — > D°K~) and A(B~ — > D°K~) amplitude ratio gives tb and 5b, and the 
values are shown in Table IVHI and compared with experiment. The tb parameter [c.f. 
Eq. (JTJ] governs the strength of the interference effect that is essential for 1^3/7 determination 
in the GLW, ADS and DK Dalitz methods. We see that our r B = 0.09 ± 0.02 prefers the 
BaBar result over the Belle result, while 5b is in agreement with both BaBar and Belle 
results. Our r# value is also in good agreement with the fit from group, obtained by 

using all three methods of GLW, ADS and DK Dalitz analysis. The smallness of tb implies 
that we need more B data to determine 03. 



It is interesting to note that the ratios y/B(D°7v°)/B(D+7i-), \ B(D°K°)/B(D+K~) and 



B(D°K~)/B(D~7T ) are enhanced by roughly 2.44, 2.85 and 1.47, respectively, from their 
factorization values. These enhancements are sometimes interpreted as 0, Q] enhancement 
in 1 02/01 1 . However, the measured ratios are nonuniversal. In particular, we would have a 
larger tb if we straightforwardly apply the 1 02/01 1 ratio from the Dtc and DK systems 15 1. 
So let us try to understand the smallness of r B or, equivalently, the smallness of the en- 



hancement in JB(D°K-)/B(D-7T°). 
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TABLE VII: Factorization anc 
the experimental results 



23, 



24 



FS I results on re, 5b with \V u b\ = 3.67 x 10 3 , and compared to 
. For the phase convention adopted see footnote |SJ 



Expt 



fac 



FSI 



r B 



>B 



0.21 ± 0.08 ± 0.03 ± 0.04 (Belle) 
< 0.19 (90% CL) (BaBar) 
0.10 ±0.04 (UT/ ft ) 
-23° ± 19° ± 11° ± 21° (Belle) 
-66° ± 41° ± 8° ± 10° (BaBar) 



0.07 ±0.03 



0.09 ± 0.02 



T(19.9i 2 Ji) 



In the factorization approach we have 

VubVcs 



fac 
' B — 



V c bV us 



eSY 2 - cff * 0.07 (19) 

with = (m 2 B - m 2 D )f K F BD (m 2 K )/[(m 2 B - m 2 K )f D F BK (m 2 D )} ~ 1.08. This r| c agrees 

with the common estimation. In Table IVHI we see that re is indeed enhanced from its 
factorization value but the enhancement is mild. Although the D°K~ rate is enhanced by 
70% from its factorization value, the amplitude is only enhanced by 30%. The rate of the 
color-allowed D°K~ mode is not affected by the quasi-elastic rescattering. Consequently, 
rs does not differ much from its factorization prediction. Note that although the predicted 
D°K~ rate has a large error, the error in re is much reduced as it is a ratio and, furthermore, 
a ratio of amplitudes. 

It is instructive to compare B~ — > D°K~ with the color suppressed modes B° — ► D°tt° 
and D°K°. Before rescattering, i.e. at factorization level, we have 



4/ eff 

D°7T ' 



0.13, 



A f D , V2afX L 



D+w - V-6«l ^DlT 

0.23, 



D°K° _ a 2 



V2 Y DK % ^ 0.46 (20) 



aL 



V- n: , , af 



D K 



At 



af 



with X DP = (m| - m 2 D )f P F BD (m 2 P )/[(m 2 B - m 2 P ) f D F BP (m 2 D )], which gives X Dn ~ 1.38, 

and Y DK = fDF^ K (m 2 D )/[f Ds F Bn (m 2 D J} ~ 1.33. The rates of the sources D+n' and 

are much larger than the corresponding color suppressed modes, hence the effects of FSI are 
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prominent. For D + ti^ 1 it can be traced to the source mode being relatively enhanced by 
Xon, while the color suppressed D°7r° is further suppressed by l/y/2 in the 7r° wavefunction. 
In the case of B~ — > D°K~, however, the situation is reversed. The color-suppressed 
mode D°K~ is relatively enhanced by Ybk-, while the D~tc° source mode receives the l/y/2 
wavefunction suppression. Thus, the factorization rate for D°K~ differs less from the D~tt° 
source rate, and the enhancement of tb through FSI is mild. Since the FSI effect is not simply 
multiplicative [see Eqs. (JT7jl and (jl8j) ]. we do not have a simplistic universal enhancement 
in \a 2 /ai\. 



C. Discussion 



Let us first offer some remarks on the fitted results of Tables IIHI and IIVI The effective 
Wilson coefficients af = 0.92+°;^ [F B D (0)/0.67] and af = 0.22ig;g [F^(m 2 D )f 0.28] ob- 
tained in our fit to DP data agree well with IoqI = 0.26 ± 0.02 from fit to B — > J/ibK 
data |42j, and with the range of af ~ 1, af ~ 0.2-0.3 from various modes [42J, |43J. The 
ratio af /af ~ 0.24 is close to the one used in our previous analysis 11]. 

We find from data that a small mixing angle |r| <C 1 is preferred. Thus, the approximation 
of treating 77 as r] 8 and ignoring r\\ taken in o ur p revious analysis [llj is basically valid, and 
the FSI phases S and 9 correspond to those in |ll| . In fact, the values for the phases S ~ ±62° 
and 9 ~ 24° are consistent with the previous results of ±48° and ±25°, respectively. We 
note that, as a consequence of the smallness of r, the phase a is less constrained from the 
present D 0r q^ data. 

From the r[ values given in Table we see that exchange rescattering is dominant over 
annihilation rescattering. Comparing f[ and f[ with r[, which should be identical in the U(3) 
limit, we see that U(3) is not a very useful limit for these modes. 

In the present work the overall sign of FSI phases cannot be determined. We may obtain 
some information comparing to other work. For example, a pole model calculation ^| 
with some inputs, such as the Wilson coefficients a\ s 2 and cut-offs for form factors in strong 
interaction vertices, or the pQCD approach can give D°h° rates in good agreement with 
D°h° data and our results. For the system, we have 

C - E\ -V2A D07T o „ 003T91 o, ; (C - E\ _ -\/2A^o ^ 



T + Cj Dn A D o„- \T + E J Dn A D+ 



°- 33eT9 ' 7^ = 7 - 0-43e^^, (21) 
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where negative (positive) phase corresponds to the case of positive (negative) 5, a and 9. 



Comparing to (C - E)/(T + C) = 0.33 e" 50 ^ |18(, 0.34e" 92Oi Q and (C - £)/(T + £) = 
0.40 e" 670i Q, 0.32 e" 111 ^ Q, the case of positive 5, a and 9 phases is preferred. It is 
interesting to note that, since the amplitudes for color suppressed modes in the Dtc system 
are fed dominantly from the same amplitude A D + W ~, by neglecting contribution from a 2 , we 
have 

C — E\ iri — ir' -rfw / C — E\ ir' t«q°,- 

5 0.41e T65 \ ~ ■ — r— — ■ — ~ 0.49e T89 \ 



T + Cj Dw l + ir' + ir' e ' \T + Ej Dn l + ir' +ir' a 

(22) 

which give estimations within 30% errors. 

In the extraction of 03 from B — > DK , D K decays the sign of the strong phase 5b can 
be determined because one is making a CP violation study. The above preferable case of 
positive FSI phases 5, o and 9 leads to a negative 5b (c.f. Table rvTT|) . which is supported 
by the measured central values. 

Let us turn to the question of r# determination, and compare with other approaches. 

We fitted six parameters, two effective Wilson coefficients, three FSI phases and one 
mixing angle, from rates of six modes, D°tt~ , D + n~ , D tt , DfK~, D°r] and D°i]'. These 
parameters are fully determined as the number of unknowns equals that of input and, con- 
sequently, the errors in parameters are propagated from data errors. To keep the above 
features we do not include more modes, such as the three DK modes as input. In principle, 
we can also include them in the fit. However, as discussed already after Eq. ([180. the tb 
obtained in the new fit should be consistent with the one given here within errors. 

One can extract the topological amplitudes from the D°h° data as in Ref. Q|. However, 
it is not clear how to apply these amplitudes to the DP system as the two topologies are 
not identical; for example, the form factor dependence in tf and Tf in Eq. (JT5j) are different. 
Furthermore, the annihilation amplitude, which turns out to be non-negligible in the DP 
system, cannot be extracted in the D°h° system. 

In pQCD approach the D°h° rates are explained through the enhancement in C from the 
incomplete cancellation in the non-factorization contribution Q|. Similar mechanism may 
lead to an enhanced Tb- Recently, a calculation in the pQCD approach gives r# = 0.093 [26], 
which is close to our result. 

Color allowed B — > Dti modes can be calculated in the framework of QCD factoriza- 
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tion [44J, but the color suppressed decay amplitudes with the emission of a heavy meson 
cannot. A process dependent a 2 approach is used j], Q. For a comparison to the present 
approach, see Ref. 11]. In this vein, soft colinear effective theory (SCET) have received 
some attention lately. Although SCET cannot predict the D°h° rate, the prediction on the 
simila rity of D*°h° and D°h° phase and the B(D°7]') j B(D°i]) ratio agree with data within 
error |l3j]. It would be interesting to see the SCET prediction on r#. We note that the 
strong phases extracted form D*°h° in our FSI approach are similar to those in D°h° jll|. 
in agreement with SCET. Consequently, r n in the D*K~ system could be similar to that 
presented here. Indeed, tb = 0.09 ± 0.04 [2£| in the D*K~ system is given by the UT/ it 
group and is close to our estimation for B — > DK~ . 



IV. CONCLUSION 

We study quasi-elastic rescattering effects in B — > DP, DK and D K modes. The 
updated B° — > D°tt , DfK~, D 0r q, D°7]' data, together with D + n~ and D°n~ , are used 
to extract a\ 2 and four rescattering parameters. We find the effective Wilson coefficients 
af ff ~ 0.92, a| ff ~ 0.22, the strong phases 5 ~ 62°, 9 = 24°, a ~ 127° and mixing angle 
t ~ 2°. The values of 5 and 9 are close to our previous ignoring Drji. The 

smallness of r implies small mixing of D°i]i with other DP modes, hence our previous 
approximation is valid. The predicted B~ — > D°K~ and B° — > D + K~, D°K° rates are in 
agreement with data. The formalism can be applied to B — > D K modes, and the rates 
for D°K~, D~K°, D~tt°, D~r) and D~r\ modes are predicted. In particular, we predict 
tb = 0.09±0.02, which agrees with the UT/ it extraction 25\ and a recent pQCD result 
Our Tb value prefers the lower value of the BaBar experiment and disfavors the Belle result, 
extracted from the ^3/7 fit to B~ {D°, D°}K~ data using the DK Dalitz method. 
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